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ABSTRACT 

We present fully-sampled maps of 461 GHz CO J = 4 ^ 3, 807 GHz GO J = 7 ^ 6, 
and 492 GHz [Gi] ^Pi —> ^Pq emission from the inner 3 degrees of the Galactic Gen- 
ter region taken with the Antarctic Submillimeter Telescope and Remote Observatory 
(AST/RO) in 2001-2002. The data cover —1?3 < £ < 2°, —0?3 < b < 0°. 2 with 0.5' 
spacing, resulting in spectra in 3 transitions at over 24,000 positions on the sky. The 
CO J = 4 —> 3 emission is found to be essentially coextensive with lower-J transitions 
of CO. The CO J = 7 —!■ 6 emission is spatially conhned to a far smaller region than 
the lower-J CO lines. The [Ci] ^Pi —> ^Pq emission has a spatial extent similar to the 
low-J CO emission, but is more diffuse. Bright CO J = 7 ^ 6 emission is detected in 
the well-known Galactic Genter clouds Sgr A and SgrB. We also detect CO J = 4—> 3 
and CO J = 7 ^ 6 absorption from spiral arms in the galactic disk at velocities near 
0 kms“^ along the line of sight to the Galactic Center. Analyzing our CO J = 7 —> 6 
and CO J = 4 —> 3 data in conjunction with J = 1 —> 0 ^^CO and ^^CO data previously 
observed with the Bell Laboratories 7-m antenna, we apply a Large Velocity Gradient 
(LVG) model to estimate the kinetic temperature and density of molecular gas in the 
inner 200 pc of the Galactic Genter region. We show maps of the derived distribution 
of gas density and kinetic temperature as a function of position and velocity for the 
entire region. Kinetic temperature was found to decrease from relatively high values 
(>70 K) at cloud edges to low values (<50 K) in the interiors. Typical gas pressures 
in the Galactic Genter gas are n(H 2 ) • Tkin ~ 10®’^ Kcm“^. We present an {£, b) map of 
molecular hydrogen column density derived from our LVG results. 

Subject headings: Galaxyxenter — Galaxy:kinematics and dynamics — ISM:atoms — 
ISM:general — ISM:molecules 
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1. Introduction 

Much has been learned about dense gas in the Galactic Center region through radio spec¬ 
troscopy. Early observations of F{2 2) OH absorption (Robinson et al. 1964; Goldstein et ah 

1964) suggested the existence of copious molecular material within 500 pc of the Galactic Center. 
This was confirmed by detection of extensive J = 1 ^ 0 ^^CO emission (Bania 1977; Liszt &: Bur¬ 
ton 1978). Subsequent CO surveys (Bitran 1987; Stark et al. 1988; Oka et al. 1998; Bitran et al. 
1997) have measured this emission with improving coverage and resolution—these surveys show a 
complex distribution of emission, which is chaotic, asymmetric, and non-planar; there are hundreds 
of clouds, shells, arcs, rings, and filaments. On scales of 100 pc to 4kpc, however, the gas is loosely 
organized around closed orbits in the rotating potential of the underlying stellar bar (Binney et al. 
1991). Some CO-emitting gas is bound into clouds and cloud complexes, and some is sheared by 
tidal forces into a molecular inter-cloud medium of a kind not seen elsewhere in the Galaxy (Stark 
et al. 1989). This diffuse inter-cloud medium appears in absorption in E(2 —> 2) OH (McGee 1970; 
Robinson & McGee 1970), in (lio ^ In) H 2 CO (Scoville, Solomon, & Thaddeus 1972), and in 
J = 0 —> 1 HCO^ and HCN (Linke, Stark & Frerking 1981). In contrast, the clouds and cloud 
complexes are dense, as they must be to survive in the galactic tide, and they appear in spectral 
lines which are tracers of high density (n(H 2 ) > 10^ cm“^), such as NH 3 (1, 1) (Glisten, Walmsley, 
& Pauls 1981) and GS J = 2 —> 1 (Bally et al. 1988). The large cloud complexes, Sgr A, SgrB, and 
Sgr C, are the among the largest molecular cloud complexes in the Galaxy (M > 10®'^ Mq). Such 
massive clouds must be sinking toward the center of the galactic gravitational well as a result of 
dynamical friction and hydrodynamic effects (Stark et al. 1991). The deposition of these massive 
lumps of gas upon the center could fuel a starburst or an eruption of the central black hole (Genzel 
& Townes 1987). 

As prelude to further study of the Galactic Genter molecular gas, we would like to determine 
its physical state—its temperature and density. This involves understanding radiative transfer in 
CO, the primary tracer of molecular gas. Also useful is an understanding of the atomic carbon 
lines, [Ci], since those lines trace the more diffuse molecular regions, where CO is destroyed by UV 
radiation but H 2 is still present. 

The J = 1 ^ 0 ^^CO line is often optically thick. Its optical depth can be estimated by 
studying its isotopomers, ^^CO and C^®0. In the Galactic Center region, ^^CO is ~ 24x less 
abundant than ^^CO (Penzias 1980; Wilson & Matteucci 1992), and C^®0 is ~ 250x less abundant 
than ^^CO (Penzias 1981). Since the radiative and collisional constants of all the isotopomers are 
similar, the ratio of optical depths in their various spectral lines should simply reflect their relative 
abundances. Where the lines are optically thin, the line brightnesses should be in the same ratio 
as the isotopic abundances; where the lines are thick, deviations from the abundance ratios are a 
measure of optical depth. Bally et al. (1987, 1988) and Stark et al. (1988) produced fully-sampled 
surveys of ^^CO and ^^CO in the Galactic Center region. They find that the ratio of the ^^CO 
J = 1 ^ 0 to ^^CO J = 1 ^ 0 line brightness temperatures is typically 10 ± 2 

in Galactic Center gas that is far from dense cloud cores. This indicates much of the Galactic 
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Center emission is only moderately thick ~ 2), especially in comparison to the galactic 

disk outside 3kpc radius, where ~ 6 (Polk et al. 1988) is smaller, even though the 

isotope ratio ~ 40 (Penzias 1980) is larger. Heiligman (1982) and Dahmen et al. (1998) 

made surveys in J = 1 ^ 0. These show J = 1 ^ 0 to J = 1 ^ 0 line 

brightness temperature ratios which vary from 40 to over 200, with typical values 

near 70, indicating values of which vary from 3 to less than 1, while the core region of SgrB2 
shows tI^q 10 . 

Determining the excitation temperature of CO works best if emission lines from several J 
levels have been measured. Lacking such observations, what is often done is to use the brightness 
temperature of the ^^CO J = 1 —> 0 line as a lower limit to the excitation temperature of the J = 1 
state, Tex,j=i. This estimate can be misleading, because the emission may not fill the telescope 
beam, diluting the brightness temperature and causing it to be many times smaller than Tex,j=i; 
as will be apparent from the data to be presented here, this is the usual case for gas in the Galactic 
Center region. 

Moving up the energy ladder, Sawada et al. (2001) surveyed the Galactic Genter region in 
^^GO J = 2 ^ 1. They compare their data to the J = 1 ^ 0 data of Bitran et al. (1997) and 
find ~ ^ 0.01, with little spatial variation. What this means is that almost all 

the CO in the Galactic Center region has low-J states which are close to local thermodynamic 
equilibrium (LTE), so that the excitation temperatures rex,j of those states are all close to the 
kinetic temperature, Tkin, and the ratio of line brightnesses for transitions between those states are 
near unity and therefore independent of Tkin (cf. Goldreich & Kwan 1974). LTE in the low-J states 
of CO does not occur under all circumstances in the interstellar medium, but it is very common 
and appears to be the rule for Galactic Center gas. Eor each value of T^in and n(H 2 ), there will, 
however, be some value of J above which all higher-J states fail to be populated, because their 
Einstein A coefficients (which increase as J^) are so large that the collision rate at that value of T^in 
and n(H 2 ) cannot maintain those states in LTE, and they must therefore be subthermally excited, 
i.e., Tex,j << Jkin- The brightness temperature of the J —> J — 1 line from those states will be 
significantly less than that of the lower-J states, and the line ratios will be much 

smaller than unity. Unlike the low-J states, the value of the line ratios from those higher-J states 
will vary from place to place, depending on Tkin, n(R 2 ), and radiative transfer effects. 

Higher still on the energy ladder, Kim et al. (2002) used the AST/RO telescope to survey a 
strip at 6 = 0° in ^^CO J = 4 ^ 3 and ^^CO J = 7 ^ 6 . They found that even the T 4 ^ 3 /T/^q ratio 
is not far from unity and shows little spatial variation. In contrast, the distribution of J = 7 ^ 6 
line emission was found to be markedly different from the lower-J transitions. Temperatures and 
densities could therefore be calculated as a function of position and velocity using the varying value 
of and an estimate of from the Bally et al. (1987, 1988) and Stark et al. (1988) 

data. In the current paper, we extend this work to a fully-sampled {£, b) map of the Galactic Genter 
region, and estimate kinetic temperature, Tkm, and density, n(H 2 ), throughout our mapped area. 
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Ojha et al. (2001) used AST/RO to make a coarse survey of the Galactic Center region in the 
(^Pi^^Pq) line of [Ci] at 492 GHz. They find that [Ci] is distributed approximately like CO, but 
that the ratio Tqi/TI‘^q is smaller in the Galactic Center region than in the disk outside 3kpc. 
This is attributable to the smaller average optical depth of the ^^CO J = 1 —> 0 in this region. The 
absence of strong spatial variations of Tqi/TI'^q across their map indicates that the amounts of 
CO-emitting and [C i]-emitting gas are approximately proportionate across the map. In the current 
paper we extend this work to a more detailed, fully-sampled (£, h) map in [Ci]. 

In § 2, we describe the observatory and the observations. In § 3, the full data set in each of 
the three observed transitions is presented in the form of spatial-spatial maps integrated over the 
full velocity range, velocity-channel maps, and spatial-velocity maps at several galactic latitudes. 
In §4, we present a Large Velocity Gradient (LVG) model of the radiative transfer in CO, which 
allows a determination of kinetic temperature, Tkin, and density, n(H 2 ). In §5, we present our 
conclusions. 


2. Observations 

The observations were performed during the austral winter seasons of 2001 and 2002 at the 
Antarctic Submillimeter Telescope and Remote Observatory (AST/RO) located at 2847 m altitude 
at the Amundsen-Scott South Pole Station. This site has very low water vapor, high atmospheric 
stability and a thin troposphere making it exceptionally good for submillimeter observations (Cham¬ 
berlin, Lane, & Stark 1997; Lane 1998). AST/RO is a 1.7m diameter, offset Gregorian telescope 
capable of observing at wavelengths between 200^um and 1.3mm (Stark et al. 1997, 2001). A 
dual-channel SIS waveguide receiver (Walker et al. 1992; Honingh et al. 1997) was used for si¬ 
multaneous 461-492 GHz and 807 GHz observations, with double-sideband noise temperatures of 
320-390K and 1050-1190K, respectively. Telescope efficiency, rji, estimated using moon scans, 
skydips, and measurements of the beam edge taper, was 81% at 461-492 GHz and 71% at 807 GHz. 
Atmosphere-corrected system temperatures ranged from 700 to 4000 K at 461-492 GHz and 9000 
to 75,000K at 807GHz. 

A multiple position-switching mode was used, with emission-free reference positions chosen at 
least 60^ from regions of interest. These reference positions were then shared by a strip of points 
at constant galactic latitude which were observed five at a time. This mapping mode caused each 
point in the map to be observed for 60 s per pass through the map. In an attempt to obtain uniform 
noise over the entire region of interest, subregions were reimaged as often as required. 

Emission from the CO J = 4 ^ 3 and CO J = 7 ^ 6 lines at 461.041 GHz and 806.652GHz, 
together with the [Ci] ^Pi ^ ^Pq and [Ci] ^P 2 ^ ^Pi lines at 492.262 GHz and 809.342 GHz, was 
imaged over the Galactic Center region —1(3 < £ < 2°, —0(3 < b < 0°2 with 0.5' spacing in £ and 
b] i.e., a spacing of a half-beamwidth or less. Smaller selected areas were also observed with longer 
integration times in the [Ci] ^Pi ^ ^Pq line. Maximum pointing errors were no larger than l', and 
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the beam sizes (FWHM) were 103-109'^ at 461-492 GHz and 58'^ at 807GHz (Stark et al. 2001). 
To facilitate comparison of the various transitions, the data were regridded onto a 0.25^ grid and 
smoothed to a FWHM spatial resolution of ‘2! with a Gaussian filter function. 

Two acousto-optical spectrometers (AOSs; Schieder, Tolls, & Winnewisser 1989) were used 
as backends. The AOSs had 1.07 MHz resolution and 0.75 GHz effective bandwidth, resulting in 
velocity resolution of 0.65kms~^ at 461 GHz and 0.37kms~^ at 807GHz. To facilitate comparison, 
the data were then smoothed to a uniform velocity resolution of Ikrns”^. The high frequency 
observations were made with the CO J = 7 —> 6 line in the lower sideband (LSB). Since the 
intermediate frequency of the AST/RO system is 1.5 GHz, the [Ci] ^P 2 ^ ^Pi line appears in 
the upper sideband (USB) and is superposed on the observed LSB spectrum. The local oscillator 
frequency was chosen so that the nominal line centers appear separated by lOOkms”^ in the double¬ 
sideband spectra. The standard chopper wheel calibration technique was employed, implemented 
at AST/RO by way of regular (every few minutes) observations of the sky and two blackbody loads 
of known temperature (Stark et al. 2001). Atmospheric transmission was monitored by regular 
skydips, and known, bright sources were observed every few hours to further check calibration and 
pointing. At periodic intervals and after tuning, the receivers were manually calibrated against a 
liquid-nitrogen-temperature load and the two blackbody loads at ambient temperature and about 
100 K. The latter process also corrects for the dark current of the AOS optical CCDs. The intensity 
calibration errors became as large as ±15% during poor weather periods. 

Once taken, the data in this survey were reduced using the COMB data reduction package. 
After elimination of scans deemed faulty for various instrumental or weather-related reasons (< 2% 
of the total dataset), linear baselines were removed from the spectra in all species by excluding 
regions where the CO J = 1 —> 0 spectra showed emission greater than > IK. This allowed 
known emission in the Galactic Center region to be readily excluded from the baseline fitting 
procedure and was generally sufficient. In a few cases, usually due to higher than average Trms for a 
given reduced spectrum, this method fails and artifacts (e.g., vertical lines in the longitude-velocity 
maps) appear. 

While the original intent was to make Trmg as uniform as possible across the entire map, this 
was not always possible. For the CO J = 4 ^ 3 transition, Ti-ms in lkms“^ wide channels with 2' 
spatial smoothing is on average < 0.3K except in the region l!8 > £ > l(5 where Trms ^ 0.8K. 
The [Ci] ^Pi ^ ^Pq transition has Ti-ms ^ 0.5 K in lkms“^ channels for the central region of 
1°0>£> -0?5, < 1.0 K for £>1°0, and < 2K for £ < -0(5. Finally, for CO J = 7 ^ 6 

(ignoring the occasional baseline feature), Trms ^ 0.8K in lkms“^ channels for £ > 0.7 and 
—0.5 > £> —0.8, and Trms ^ 2K elsewhere. 
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3. Data Presentation 

Sample spectra at the respective positions of peak CO J = 7 ^ 6 emission toward the Sgr A 
{£= 0?00, b= -0?07), SgrB {£= 0?66, b= -0?05), SgrC {£= -0°45, b= -0°20) (CO J = 4 ^ 3 
emission peak), and £ ~ 1°3 {£= 1°25, b= —0?05) molecular complexes are shown in Fig. 1. The 
CO J = 4 ^ 3, CO J = 7 —> 6, and [Ci] ^Pi ^ ^Pq data are from the AST/RO survey, while the 
CO J = 1 ^ 0, ^^CO J = 1 ^ 0, and CS J = 2 ^ 1 data are from the Bell Laboratories (BL) 
7-m telescope (Stark et al. 1988; Bally et al. 1987, 1988). All of the spectra are from datacubes 
smoothed to 2’ resolution. 

The CO J = 4 ^ 3 profile toward Sgr B resembles the CO J = 1 ^ 0 profile, with similar 
linewidth and a prominent self-absorption feature at ulsr = 60 kms“^. The total velocity extent is 
somewhat smaller in CO J = 4 ^ 3 on the negative velocity side. In contrast, the [Ci] ^Pi ^ ^Pq 
and CO J = 7 ^ 6 lines, as well as the ^^CO, show peak emission at the self-absorption velocity, 
suggesting these lines are less optically thick. It is important to note that the strong feature at 
negative velocities in the CO J = 7 ^ 6 spectrum toward Sgr B is due to superposed emission in 
the 809 GHz (^P 2 ^^Pi) line of [Ci] in the image sideband. 

Fig. 2 presents spatial-spatial {£, b) maps integrated over velocity for the three transitions 
observed with AST/RO^ and, for comparison, the three transitions observed at the BL 7-m. All 
six maps have been smoothed to the same 2' spatial resolution. The regions with £ > 0?9 and 
£ < —0?5 in the ^^CO J = 1 ^ 0 map were observed with sparser sampling than the rest of the 
map (Bally et al. 1987). The most striking result is that CO J = 4 ^ 3 emission in the Galactic 
Genter region is essentially coextensive with the emission from the lower J transitions of GO. This 
contrasts sharply with the outer Galaxy where GO J = 4 ^ 3 emission is rather less extensive 
than GO J = 1 —> 0. In all six maps, the brightest emission occurs primarily at negative latitudes. 
Four major cloud complexes are seen in the GO J = 1 ^ 0, GO J = 4 ^ 3, and ^^GO J = 1 ^ 0 
maps, from left to right: the complex at £ ~ 1°3, the SgrB complex near £ ~ 0°7, the Sgr A cloud 
near 7 ~ 0°0, and the SgrC cloud near {£ ~ —0°45, b ~ —0°2). As noted by Kim et al. (2002), 
the CO J = 7 —> 6 emission is much more spatially confined than the lower-J CO transitions. 
In contrast, the [C i] emission is comparable in spatial extent to the low-J CO emission, but its 
distribution appears somewhat more diffuse (less peaked). The SgrC cloud is much less prominent 
in the [Cl] map than in the other five transitions. The noise in the [Ci] map at £ < —0)’7 is greater 
than that in the rest of the map due to shorter integration times. 

Figs. 3-5 present spatial-spatial maps in successive lOkms”^ wide velocity ranges in the 
CO J = 4 ^ 3, [Ci] ^Pi ^ ^Pq, and CO J = 7 ^ 6 lines, respectively. For each spectral line, 
velocity channels centered at velocities from —165kms~^ to 4-185 kms~^ (see upper left corner of 


^Complete FITS datacubes of the three AST/RO datasets are archived with the electronic edi¬ 
tion of this paper. As new data becomes available updated versions will be posted at http://cfa- 
www .harvard .edu/~adair/A S T_RO / ab c. ht ml. 
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each map) are displayed. In each map, the integrated intensity has been divided by 10 km s ^ so that 
the intensity color scale is a reasonably good indicator of corrected average antenna temperature, 
(^ 1 ). 

Spatial-velocity maps at 12 different values of galactic latitude are shown in Figs. 6-8 for 
CO J = 4 —> 3, [Ci] ^Pi —> ^Pq, and CO J = 7 ^ 6, respectively. All major structures in the 
galactic core can be identified in these maps, including Sgr A, Sgr B, Sgr C, and the £ ~ l)3 complex. 
The 300 pc molecular ring shows up prominently in both the [C i] and CO J = 4 ^ 3 {i- v) maps, 
running in a straight, continuous line from {i ~ 2°,ulsr ~ 220 kms“^) to {£ ~ —0°9,ulsr ~ 
120kms“^), and from {£ ~ 1°5,ulsr ~ —20kms“^) to {£ ~ —0°9,ulsr ~ —180kms“^). There is 
little or no evidence that the CO J = 7 ^ 6 line is excited into emission in the 300 pc ring. Also 
absent in the CO J = 7 ^ 6 line is any trace of the 3 kpc arm, which is prominent in CO J = 1 ^ 0 
emission. The 3 kpc arm must contain little dense, warm gas. Foreground absorption by spiral arms 
in the galactic disk is seen at velocities near Okrns”^ in CO J = 4 ^ 3 and [Ci] ^Pi ^ ^Pq. A 
hint of this absorption feature is seen in CO J = 7 ^ 6. 

4. LVG Model 

The large velocity gradient (LVG) approximation (Goldreich & Kwan 1974) simplifies radiative 
transfer analysis of molecular lines. Imagine a small volume within a molecular cloud, in which a 
molecule emits a photon. The volume is sufficiently small that all velocities within it are thermal, 
and temperature and density are constant. It is surrounded by other small volumes which have 
similar temperature and density, and whose internal velocities are also thermal, but the velocities 
of these other volumes are different from that of the volume which emitted the photon because of 
velocity gradients within the molecular cloud. In LVG, it is assumed that if the photon is going 
to be absorbed, rather than escaping from the cloud, it can only be absorbed in nearby volumes, 
those volumes whose velocities are close to the velocity of the emitting volume. In this case ‘close’ 
means close enough that the thermal linewidths overlap. As the emitted photon travels away from 
its point of origin, it passes through nearby volumes where it has a chance to be absorbed because 
the thermal motion of some of the molecules causes them to have the same velocity as the emitting 
molecule. At some distance from the point of emission, this will no longer be true because of 
velocity gradients, and in LVG it is assumed that the photon then escapes from the cloud. The 
velocity gradient is ‘large’ in the sense that the region of possible absorption is small compared to 
the cloud as a whole, so it can be assumed that in effect the photon was absorbed at its point of 
origin. This achieves a great simplihcation in making the radiative transfer problem entirely ‘local’, 
while allowing for the possibility of absorption. The probability that a photon will escape from 
the cloud and become observable then depends only on the local physical and chemical properties 
of the emitting gas and on the value of the velocity gradient. The emitting cloud can then be 
modeled using only a few parameters, yet the radiative transfer is realistic enough to compare to 
observations. The parameters of the model can then be adjusted to fit the data, and an estimate 



of those parameters is thereby achieved. Ossenkopf (1997) has shown that such estimates are 
robust, in the sense that the parameter values derived are often approximately correct, even in 
circumstances where the assumptions of the LVG approximation are violated. 

We will use this LVG methodology to estimate the kinetic temperature, Tkm, and the number 
density of molecular hydrogen, n(H 2 ), throughout the Galactic Genter region. Due to the high 
velocity dispersions characteristic of the Galactic Genter, the LVG approximation is most likely 
valid over much of the mapped region. The LVG approximation does not apply to some of our 
data: the foreground absorption by spiral arms is clearly not local to the emitting gas in the 
Sgr A cloud. It should be kept in mind that the LVG approximation is an Ansatz which allows 
us to estimate the physical properties of the Galactic Genter gas in what otherwise would be an 
untenably complex modeling problem. 

Our cloud model was developed by M. Yan and S. Kim. It has plane-parallel cloud geometry. 
It uses GO collisional rates determined by Turner (1995) and uses newly-derived values for the H 2 
ortho-to-para ratio (~ 2) (Rodriguez-Fernandez et al. 2000) and for the collisional quenching rate 
of CO by H 2 impact (Balakrishnan, Yan, & Dalgarno 2002). The model has two input parameters: 
the ratio of ^^CO to ^^CO abundance, and the ratio Y(C'0)/W, where X{CO) is the fractional 
CO abundance parameter and W is the velocity gradient. As discussed in §1, the abundance ratio 
^^CO/^^CO is 24 in the Galactic Center region (Penzias 1980; Langer & Penzias 1990; Wilson & 
Matteucci 1992; Langer & Penzias 1993). We will use a value X{CO)/W = 10“^'® pckm“^s, 
assuming that the ^^CO/H 2 ratio is 10“^ and the velocity gradient within the Galactic Center 
gas is a uniform 3 kms~^pc“^. Dahmen et al. (1998) estimated that W is 3 kms~^pc“^ to 
6 kms“^pc“^, and indeed these are typical slopes of position-velocity features in Galactic Center 
maps. There is, however, no reason to suppose that a single value of W applies throughout the 
Galactic Center region, or even that W is constant within a single cloud. This is a weak point 
in the analysis, because in the LVG analysis n(H 2 ) oc (W)^'®, when all other parameters are held 
fixed. 

For each observed point, we take the brightness temperature ratios and T/^^q/T/^q, 

using the same methodology as Kim et al. (2002), to determine Tkm and n(H 2 ). Fig. 9 is a smoothed 
representation of the relationships generated by our LVG model. Note that our model, using the 
ratio, is particularly sensitive to variations in density and temperature near ~ 40K 
and n(H 2 ) ~ 10^ cm“^, and that its range of validity extends a factor of ~ 20 in density and ~ 4 in 
temperature around these values. Much of the Galactic Center gas falls within this range; the gas 
which does not could be studied using various other transitions of the CO isotopomers, or using 
other species such as [Ci] and CS. 

Fig. 10 shows the output of our LVG analysis for Tkin as spatial-velocity maps for 6 different 
values of galactic latitude, while Fig. 11 does the same for n(H 2 ). White regions on the maps 
indicate areas where either the LVG inversion did not converge or where there is missing data from 
the maps of the four transitions used. The LVG model strains to fit in regions where 
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so small it approaches the value of the isotope ratio (1/24). Where < 0.1, that is, where 

the J = 1 ^ 0 CO line is less optically thick, the model assigns high kinetic temperatures (cf. the 
bottom of Fig. 9). Tkin starts high in the lower density regions. This is physically reasonable since 
the edges of these complexes are likely to be of lower density and externally heated. Moving into the 
interior of the complexes, we find a surprisingly uniform temperature of T^in ~ 50 K and density of 
n(H 2 ) ~ 10^'® cm“^. This indicates a typical pressure n(H 2 ) • Tkin ~ 10^'^ Kcm“^. Exceptional are 
the density peak at SgrB and the foreground material along ulsr — Okrns”^. The breakdown of 
LVG assumptions in the foreground absorbing material renders the Tkin and n(H 2 ) results suspect 
if not outright invalid near the absorbing material. 

Fig. 12 shows the estimated column density as determined by integrating the LVG spatial 
density and dividing by W. The region —60 < ulsr < 20kms“^ is excluded in order to avoid 
contamination by the foreground material for which the LVG analysis is invalid; the actual total 
column density should therefore be somewhat larger than the value shown here. Also, there are a 
few places where the LVG model fails to converge even though the lines are strong and the density 
is presumably high—the density in those spots does not contribute to our estimate of the column 
density. There is an overall factor of order unity uncertainty in the scaling of this map due to likely 
errors in our assumed value of W. The map clearly shows the density enhancement toward Sgr A 
and Sgr B and the relatively high density of the material surrounding and connecting them. In the 
closed-orbit paradigm of Binney et al. (1991), this material is in X 2 orbits. 


5. Conclusions 

1. We have mapped the inner 3° of the Galaxy in 461 GHz GO J = 4 —> 3, 807 GHz GO J = 
7^6, and 492 GHz [Gi] ^Pi ^ ^Pq emission. Sgr A, SgrB, SgrG, and the 300 pc molecular 
ring are easily identified. 

2. The GO J = 4 ^ 3 emission is found to be essentially coextensive with lower-J transitions 
of GO. The GO J = 7 ^ 6 emission is spatially confined to a far smaller region than the 
lower-J GO lines. The [Gi] ^Pi ^ ^Pq emission has a spatial extent similar to the low-J GO 
emission, but is more diffuse. 

3. Gonsistent with Kim et al. (2002), we find the line ratio is approximately constant 

and not far from unity over much of the mapped region. In contrast, the line 

ratio is found to vary significantly. 

4. For each observed point, line ratios, together with ratios, were 

used to estimate kinetic temperatures and molecular hydrogen volume densities. Kinetic 
temperature was found to decrease from relatively high values (>70 K) at cloud edges to low 
values (<50 K) in the interiors. Molecular hydrogen densities, n(H 2 ), ranged up to the limit 
of our ability to determine via our LVG analysis, ~ 10^-^ cm-3. 
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5. Typical gas pressures in the Galactic Center gas are n(H 2 ) ~ 10®'^ Kcm“^, while typical 

virial pressures are n(H 2 ) -Tvii-iai ~ 10®'® K cm“®. These values can be compared to the typical 
gas pressures in molecular clouds near the Sun ~ 10®'^Kcm“®, the typical virial pressure in 
molecular clouds near the Sun ~ 10®Kcm“®, and the ambient pressure of the interstellar 
medium near the Sun lO^Kcm ® (Dickey & Lockman 1990). 
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Fig. 1.— Spectra toward the respective positions of peak CO J = 7 —> 6 emission in the Sgr A 
{£= 0?00, b= -0?07), SgrB {£= 0?66, b= -0?05), SgrC {£= -0°45, b= -0°20) (CO J = 4 ^ 3 
emission peak)', and £ ~ 1?3 {£= 1°25, b= —0°05) clouds (as indicated at lower right in each 
frame) in 6 different transitions, as indicated by the color identifications at upper left. The 461 GHz 
CO J = 4 ^ 3, 807 GHz CO J = 7 ^ 6, and 492 GHz [C i] ^ 3p^ AST/RO 

survey {this paper), and the 115 GHz GO J = 1 —0, 110 GHz ^^CO J = 1 ^ 0, and 98 GHz 
GS J = 2 —> 1 data are from the Bell Laboratories 7-m telescope (Stark et al. 1988; Bally et al. 
1987, 1988). 
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Fig. 2.— Spatial-spatial {i, b) integrated intensity maps for the 3 transitions observed with 
AST/RO {top 3 panels) and, for comparison, the 3 transitions observed at the BL 7-m (Stark et al. 
1988; Bally et al. 1987, 1988) {bottom 3 panels). Transitions are identified at left on each panel. The 
emission is integrated over all velocities where data are available. These values of {vmin, Vmax) are: 
[Ci], (-90,150); CO(7~6), (-30,120); CO(4-3), (-150,150); CO(l-O), (-150,150); ^300(1-0), 
(—150,150); CS(2-1), (—150,150). All 6 maps have been smoothed to the same 2' resolution. This 
figure has been corrected as noted in the erratum. 
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L.1 L a.1 □ -n.i -L 


Fig. 3.— False color velocity-channel maps for the CO J = 4 —> 3 transition toward the Galactic 
Center. Each subpanel shows a lOkms”^ velocity bin, centered at the velocity shown in the upper 
lefthand corner. The integrated emission in K kms~^ has been divided by the lOkms”^ width of 
the bin so that the color scale at right indicates (T^). This figure has been corrected as noted in 
the erratum. 
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Fig. 3.— Continued 
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Fig. 4.— False color velocity-channel maps for the [Ci] ^Pi ^ ^Pq transition toward the Galactic 
Center. Each subpanel shows a lOkms”^ velocity bin, centered at the velocity shown in the upper 
lefthand corner. The integrated emission in K kms~^ has been divided by the lOkms”^ width of 
the bin so that the color scale at right indicates {T^). This figure has been corrected as noted in 
the erratum. 
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Fig. 4.— Continued 
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Fig. 4.— Continued 
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Fig. 5.— False color velocity-channel maps for the CO J = 7 —> 6 transition toward the Galactic 
Center. Each subpanel shows a lOkms”^ velocity bin, centered at the velocity shown in the upper 
lefthand corner. The integrated emission in K kms~^ has been divided by the lOkms”^ width of 
the bin so that the color scale at right indicates (T^). This figure has been corrected as noted in 
the erratum. 
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Fig. 5.— Continued 
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Fig. 6.— False color longitude-velocity maps of CO J = 4 —> 3 emission toward the Galactic 
Center. Each of the 12 panels displays emission at a different value of galactic latitude, indicated 
in the lower left corner of each panel. 
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Fig. 7.— False color longitude-velocity maps of [Ci] ^Pi ^ ^Pq emission toward the Galactic 
Center. Each of the 12 panels displays emission at a different value of galactic latitude, indicated 
in the lower left corner of each panel. 
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Fig. 8.— False color longitude-velocity maps of CO J = 7 —> 6 emission toward the Galactic 
Center. Each of the 12 panels displays emission at a different value of galactic latitude, indicated 
in the lower left corner of each panel. 
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Fig. 9.— Approximate representation of the relation between the line ratios and T^in (blue curves, 
units are K) and n(H 2 ) (red curves, units are log[n(H2)/1.0cm“^]) generated by our LVG model, 
which uses an abundance ratio ^^CO/^^CO = 24 and X{CO)/W = 10“^'® pckm“^ s. Instabilities 
in the generation of these functional relationships have been smoothed by hand. 
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Fig. 10.— False color longitude-velocity maps of Tkm as determined by the LVG model described 
in the text. Regions in white indicate areas where either spectral line data are not available or the 
LVG model did not converge. Each of the 6 panels displays Tkin at a different value of galactic 
latitude, indicated in the lower left corner of each panel. 
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Fig. 10.— Continued 
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Fig. 11.— False color longitude-velocity maps of logn(H 2 ) as determined by the LVG model 
described in the text. Regions in white indicate areas where either spectral line data are not 
available or the LVG model did not converge. Each of the 6 panels displays n(H 2 ) at a different 
value of galactic latitude, indicated in the lower left corner of each panel. 
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Fig. 11.— Continued 
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Fig. 12.— False color velocity-channel map of log[/ n{R 2 )dv/W\ as determined by the LVG 
model described in the text. n(H 2 ) is integrated over the ranges —150 < ulsr < —60kms“^ and 
20 < i'lsr < 150kms~^ in order to avoid contamination by the foreground material for which the 
LVG analysis is invalid. This value is then divided by W in order to make a map comparable to 
the expected column density in units of cm“^. 





